A note on quantum error correction with continuous variables 
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We demonstrate that continuous-variable quantum error correction based on Gaussian ancilla 
states and Gaussian operations (for encoding, syndrome extraction, and recovery) can be very useful 
to suppress the effect of non-Gaussian error channels. For a certain class of stochastic error models, 
reminiscent of those typically considered in the qubit case, quantum error correction codes designed 
for single-channel errors may enhance the transfer fidelities even when errors occur in every channel 
employed for transmitting the encoded state. In fact, in this case, the error-correcting capability of 
the continuous-variable scheme turns out to be higher than that of its discrete- variable analogues. 

PACS numbers: 03.67.Hk, 03.67.Mn, 42.50.Pq 



INTRODUCTION 



Quantum error correction is an essential tool in or- 
der to protect fragile quantum states whenever they are 
subject to errors in a quantum computation protocol, in- 
cluding unwanted decoherence effects due to interactions 
with the environment In the case of systems de- 
scribed by discrete quantum variables such as qubits, the 
theory of quantum error correction codes (QECC) is very 
advanced. Various codes have been proposed to correct 
arbitrary types of errors. In these codes, a signal state is 
encoded into a larger Hilbert space, typically including 
additional ancilla systems, resulting in an encoded, mul- 
tipartite entangled state. Local errors on some subsys- 
tems can then be corrected, as the quantum information 
remains undisturbed in the global structure of the total 
encoded state. The simplest qubit QECC are designed 
such that only single-qubit errors can be corrected, i.e., 
errors occurring in a single channel, assuming that the 
subsystems of the encoded state are sent through indi- 
vidual, independent channels. Provided the error proba- 
bilities of such single-channel errors are below a certain 
threshold (and hence, multiple-channel errors are highly 
unlikely), a protocol based upon QECC leads to a better 
performance than an unencoded scheme. 

In the case of systems described by continuous quan- 
tum variables, the situation is not so clear. Although 
there are proposals of QECC in the regime of continu- 
ous variables [2|, y, l|, the applicability of these codes, 
and to what extent they are useful, has not been fully 
understood. For example, encoding an arbitrary signal 
state into a nine-mode wavepacket code would enable 
one to correct an arbitrary single-mode error, including 
Gaussian or non- Gaussian errors (where a Gaussian error 
maps a Gaussian state back to a Gaussian state). The 
whole protocol would nonetheless rely upon only Gaus- 
sian ancilla states (squeezed states) and Gaussian op- 
erations (beam splitter transformations, homodyne de- 
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tection, and feedforward; online squeezing transforma- 
tions arc not needed Q). This would lead to very effi- 
cient implementations of QECC. However, similar to the 
qubit case, do these codes also improve the performance 
of quantum information protocols in the presence of in- 
dependent, multiple- channel errors? Or do we have to 
strictly (and rather artificially) assume that only a single 
mode in a single channel is affected by an error, while all 
the other channels are ideal? The answer to this question 
will certainly depend on the error (or noise) model taken 
into account. In fact, quite intuitively, the transmission 
of an optical mode in a Gaussian state subject to pho- 
ton loss cannot be enhanced using the Gaussian QECC, 
as the losses would occur in every channel of the en- 
coded state and the corresponding, realistic noise model 
would lack the stochastic nature of the most common 
qubit channels. More generally, it has been proven that 
Gaussian QECC are not useful to protect Gaussian states 
against Gaussian errors [f| @ . 

Here, we will focus on a simple Gaussian three-mode 
repetition code which can correct arbitrary "a;-errors" , 
i.e., x-displacements and any other errors decompos- 
able into x-displacements (including non-Gaussian er- 
rors). A code for correcting arbitrary errors includ- 
ing non-commuting x and p-errors is obtainable, for in- 
stance, by concatenating the three-mode code into a 
nine-mode code [J]. We find that for a certain class 
of (non-Gaussian) error models, the deterministic three- 
mode QECC protocol enables one to achieve a significant 
improvement of fidelity compared to the direct transmis- 
sion of the signal state. Such error models may describe, 
for instance, free-space channels with atmospheric fluc- 
tuations causing beam jitter, as considered recently for 
various non-deterministic distillation protocols [1, 0, [l(| 
(see also fill]). 



II. STOCHASTIC ERROR MODELS 

Let us consider the following error model. The in- 
put state described by the Wigner function W lrL is trans- 
formed into a new state Horror with probability 7; it re- 
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mains unchanged with probability 1 — 7. Thus, we have 

Wont{x,p) = (1 - j)W in (x,p) + jW eiror (x,p) . (1) 

In general, the Wigner functions W- m and Mirror may de- 
scribe arbitrary quantum states. Note that even in the 



case of two Gaussian states, W in and also W r , 



the 



resulting state W out is no longer Gaussian. Thus, this 
channel model describes a certain, simple form of non- 
Gaussian errors. It is a generalization of the "erasure" 
channel considered in Ref. 0] , where a coherent-state in- 
put is displaced to a vacuum state with an error probabil- 
ity of 7; otherwise the coherent state leaves the channel 
untouched. 

Similarly, as an example, we will now consider a 
coherent-state input, = |a;i + ipi), described by the 
Wigner function, 

W in (x uPl ) = -eM-2{x 1 -x 1 ) 2 -2(p 1 -p 1 ) 2 }. (2) 

Moreover, for simplicity, we assume that the effect of the 
error is just an ^-displacement by x 2 such that 



terror = Wi n (xi - X 2 ,p\) 



(3) 



Note that more general errors, including non-Gaussian 
x-errors, could be considered as well. The sign of the 
displacement error shall be fixed and known, e.g., without 
loss of generality, x 2 > 0. 



III. ENCODING AND TRANSMISSION 



Note that this encoding procedure does not require any 
online implementations of continuous-variable CNOT 
gates 0, El ; combining the signal mode with the offline 
squeezed ancilla modes at a sequence of beam splitters is 
sufficient 

Now we send the three modes through individual chan- 
nels where each channel acts independently upon every 
mode as described by Eq. |T]) with Mirror corresponding 
to an ^-displacement by x 2 . As a result, the three noisy 
channels will turn the encoded state into the following 
three-mode state, 



(7) 



Wi nc (ai,a 2 ,a 3 ) 

= (1 - 7) 3 W nc(ai, «2, 03) 

+7(1 - 7) 2 Wo„ c (xi - x 2 + ipi, a 2 , o 3 ) 

+7(1 - ~f) 2 W cnc (a 1 ,x 2 ~x 2 + ip 2 ,a 3 ) 

+7(1 - j) 2 W cnc (a 1 ,a 2 ,x 3 -x 2 + ip 3 ) 

+7 2 (1 - 7)W cnc (a; 1 - x 2 + ipi,x 2 - x 2 + ip 2 ,a 3 ) 

+7 2 (1 -7)W enc (a;i - x 2 + ip l7 a 2 ,x 3 - x 2 + ip 3 ) 

+7 2 (1 - 7)Wenc(Q!l 5 »2 ~ X 2 + ip 2 ,X 3 ~ X 2 + ip 3 ) 

+l' 3 W enc {xi -x 2 + ipi,x 2 -x 2 + ip 2l x 3 - x 2 + ip 3 ) 



Note that we assumed the same ^-displacements in every 
channel. Let us now consider the decoding procedure and 
how to extract the error syndromes by using homodyne 
detections on the ancilla modes. 



Now in order to encode the input state, we use two 
ancilla modes, each in a single-mode x-squeezed vacuum 
state, represented by 



W anc (x k ,Pk) = -exp[-2e 

7T 



+ 2 ^ 2 - 2e" 2 > 2 ] , 



(4) 



with squeezing parameter r and k = 2,3. The total three- 
mode state before encoding is 

W{ai, 0-2,0:3) = W m (xi,pi)W anc (x 2 ,p 2 )W a , nc (x 3 ,p 3 ),(J>) 



with a-j 



i-Pji 3 = 1,2,3. The encoding may 



be achieved by applying a "fritter", i.e., a sequence of 
two beam splitters with transmittances 1 : 2 and 1:1. 
The total, encoded state will be an entangled three-mode 
Gaussian state with Wigner function, 



IV. DECODING AND SYNDROME 
EXTRACTION 



The decoding procedure now simply means inverting 
the fritter, which results in 



Wdec(ai, 02,03) (8) 

= (1 - 7) 3 Wi n {xi,Pl)W anc {x 2i p 2 ) Wane (x 3 ,p 3 ) 

1 



W e nc(0l,0 2 ,a 3 ) 



(6) 
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+7(1 - 7) Wi n ( Xi ~ -j=X 2i pi 



xW anc [ Xi - \j ~X 2 ,p 2 I W anc (x3,p 3 ) 



+7(1 - 7) 2 W„ Xl 7=X 2 ,Pl 
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xW ;m! j + -^=X 2 ,p 2 ^j Wanc(^3 - ^=X 2 ,p 3 ) 
1 



+7(1 - 7) W n ( Xl - -j=X 2 ,pi 

\'W,„, ( r, + -^=X 2 ,p 2 ^j W anc (x 3 + "^* 2 '^ 3 ) 



+7 (1 - 7)Win ( Xl - -j=Xl,pi 

\ tr all! j .co - -^=X2,P2^ W anc (x 3 - -^=X 2 ,p 3 ) 



- j)W in ( Xl - ~^ X 2,Pl 
"- W'.un { -'2 - -^=X 2 ,p 2 ^ VFanc^a + -^=X 2 ,p 3 ) 



+7 (1 - 7)W / in ( Xl - — 3J 2 ,Pl 



><Wa nc ]^X 2 + \j ^X 2 ,p 2 j W anc (x 3 ,p 3 ) 

+j 3 W in (xi - V3x 2 ,pi 

xWanc (X 2 ,p 2 ) W anc (x 3i p 3 ). 

By looking at this state, we can easily see that x- 
homodyne detections of the ancilla modes 2 and 3 
(the syndrome measurements) will almost unambigu- 
ously identify in which channel a displacement error oc- 
curred and how many modes were subject to a displace- 
ment error. The only ambiguity comes from the case 
of an error occurring in every channel at the same time 
(with probability 7 3 ), which is indistinguishable from the 
case where no error at all happens. In both cases, the two 
ancilla modes are transformed via decoding back into the 
two initial single-mode squeezed vacuum states. All the 
other cases, however, can be identified, provided the ini- 
tial squeezing r is sufficiently large such that the displace- 
ments rx x 2 , originating from the errors, can be resolved 
in the ancilla states. This will always be possible when 
r — * co and/or x 2 1. Note that even without squeez- 
ing, r = 0, all those cases can be distinguished, provided 
that x 2 ^> 1 holds. However, in this case, the recovery 
displacement operations will result in larger excess noises. 
Perfect recovery with unit fidelity requires infinite squeez- 
ing. Similarly, small displacements x 2 would require suf- 
ficiently large squeezing to be resolved, e _2 ' /4 < x 2 . 

The recovery operation, i.e., the final phase-space dis- 
placement of mode 1 depends on the syndrome mea- 
surement results for modes 2 and 3 which are consistent 
with either undisplaced squeezed vacuum states ('0') or 
squeezed vacua displaced in either '+' or '— ' x-direction. 
The syndrome results for modes 2 and 3 corresponding 
to the eight possibilities for the errors occurring in the 
three channels are (0,0) for no error at all, (+,0) for an 
error in channel 1, (— ,+) for an error in channel 2, (— ,— ) 
for an error in channel 3, (+,+) for errors in channels 1 
and 2, (+,— ) for errors in channels 1 and 3, (—,0) for 
errors in channels 2 and 3, and, again, (0,0) for errors 
occurring in all three channels. 

The (unnormalized) conditional states of mode 1 de- 
pending on the syndrome measurement results x 2 and 
x 3l including suitable feedforward operations (i.e., dis- 
placements of mode 1 using the measured results), can 



be obtained by integrating Eq. {8]| over p 2 and 
Wcoud (ai| £2,^3) 
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+7 2 (l-7) 
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= -2e +2r (x 2 + V2/32 2 ) 
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We observe that in almost all cases, the feedforward oper- 
ations turn mode 1 back into the initial state up to some 
Gaussian-distributed excess noise depending on the de- 
gree of squeezing used for the encoding. The only case 
for which no correction occurs is when errors appear in 
every channel at the same time, at a probability of 7 3 . 
In this case, the initial state remains uncorrected, with 
an ^-displacement error of \/3x2- 

In the limit of infinite squeezing, r — > oo, the Gaussian 
distribution functions in Eq. © become delta functions. 
As a result, the ensemble output state of mode 1 upon 
averaging over all syndrome measurement results X2 and 
X3 (by integrating over X2 and £3) becomes 

(1 - j 3 )W iTl ( Xl ,px) + 7 3 W in (xt - s/3x 2 ,Pi) ■ (10) 

We see that a fidelity of at least 1 — 7 3 can be achieved 
(assuming 22 1; for small X 2, the fidelity would ex- 
ceed 1 — 7 3 , but those smaller X2 may be too hard to 
detect at the syndrome extraction, depending on the de- 
gree of squeezing, see below). This result implies that 
the encoded scheme performs better than the unencoded 
scheme (direct transmission with indirect = 1 — 7) for any 
< 7 < 1. In other words, by employing the quantum 
error correction protocol, the error probability can be re- 
duced from 7 to 7 3 . The continuous- variable scheme, in 
this model, is more efficient than the analogous qubit rep- 
etition code, and it does not require error probabilities 
7 < 1/2 as for the case of qubit bit-flip errors [l|. 

Consider now the regime e _2r /4 < X2 < 1/4, corre- 
sponding to small displacements below the vacuum limit. 
The resulting displacements can only be resolved pro- 
vided the squeezing is large enough. In the limit of in- 
finite squeezing r — > 00, arbitrarily small shifts can be 
detected and perfectly corrected (with zero excess noise 
in the output states corresponding to unit fidelity). In 
the regime 12 > 1, corresponding to large shifts, even 
zero squeezing in the ancilla modes (i.e., vacuum ancilla 
states) is sufficient for error identification. For r = and 
X2 3> 1, the syndrome measurements still provide enough 
information on the location of the error and, to some ex- 
tent, on the size of the error. However, for this case of 
"classical error correction" , the recovery displacements 
lead to finite excess noises in the output state after error 
correction, originating from the vacuum ancilla states. 

V. CONCLUSIONS 

Using the simple example of a three-wavepacket repe- 
tition code, we demonstrated that for certain stochastic 



error models, the continuous- variable, Gaussian proto- 
col (based on offline squeezing, beam splitter transfor- 
mations, and homodyne detection) leads to a significant 
improvement of fidelity even when the errors occur in 
every channel. In this case, the errors correspond to er- 
rors in one variable, e.g. x-displacements or any errors 
decomposable into x-displacements. The appropriate er- 
ror model is reminiscent of the most typical qubit chan- 
nels such as a bit-flip channel. In the continuous- variable 
regime, these types of stochastic errors map a Gaussian 
signal state into a non-Gaussian state represented by a 
discrete, incoherent mixture of the input state with a 
Gaussian (or a non-Gaussian) state; thus, circumventing 
the recent nogo result on Gaussian QECC ||. 

It turns out that the fidelity gain through encod- 
ing compared to direct transmission without encoding 
is larger than that for the analogous qubit scheme. In 
fact, for the three-qubit bit-flip repetition code, the 
eight-dimensional, physical Hilbert space can be divided 
into only four orthogonal, logical qubit subspaces, cor- 
responding to the four cases of no error at all and a 
bit flip on any one of the three qubits (corresponding to 
two classical syndrome bits). In the continuous- variable 
case, even those events with errors occurring on two of 
the three modes simultaneously can be unambiguously 
identified and corrected, since more error subspaces are 
available and a correspondingly larger amount of syn- 
drome information. Possible applications and extensions 
of the scheme considered here are stochastic error mod- 
els with quadratic or even cubic or higher-order ^-errors, 
gain optimizations for the recovery displacements with 
finite squeezing, higher-level repetition codes in one vari- 
able, and non-commuting ("truly quantum") errors in 
both x and p, requiring more complex codes than just 
three modes. 
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